Comment on "Glassy Transition in a Disor- 
dered Model for the RNA Secondary Structure" 

In a recent, very interesting paper, Pagnani, Parisi and 
Ricci-Tersenghi have studied the low-temperature be- 
havior of a model for RNA secondary structure. They 
claim that the model exhibits a breaking of the replica 
symmetry, since the width of the distribution P{q) of 
overlaps may converge to a finite value at T = 0. The 
authors used an exact enumeration method to obtain all 
ground states for a given RNA sequence. Because of 
the exponential growing degeneracy, only sequences up 
to length L = 256 could be studied. 

Here it is shown that, in contrast to the previous re- 
sults, by going to much larger sizes as L = 2000 the 
variance tT^(q) oc L^" '"'. This means that P{q) becomes 
a delta function in the thermodynamical limit at T = 0. 
The method used here combines the ideas presented in 
and The method is faster than the algorithm 
of 1^ because no floating-point arithmetic is necessary. 
Furthermore, the algorithm of ||^ is not exact, although 
usually true ground states are obtained. The technique 
of ^ guarantees exact ground states but is restricted to 
small sizes. 

Here, a finite number of exact ground states is selected 
randomly from the set of all ground states which is rep- 
resented by a graph. Similar to an ordinary Monte-Carlo 
simulation it has to be guaranteed that each ground state 
appears with the proper weight, i.e. with the same prob- 
ability, since all ground states have exactly the same en- 
ergy. This is ensured by the following technique: Let 
now Gij denote the set of ground states for the sequence 
[vi, . . . ,rj]. Similar to the representation of the parti- 
tion function applied in [^,0, Gij can be expressed in 
terms of ground states for smaller sequences: a ground 
state for the sequence [r^, . . . , r^] can either be a ground- 
state of [r^+i, . . . , 7'j] (if the energy is low enough), or it 
is a combination of a pair (i, fc) (fc £ (i -I- 1, . . . , j)) with 
an arbitrary ground state of [r^, . . . ,rk-i] and an arbi- 
trary ground state of [rk+i, ■ ■ ■ ,rj] (if the energy is low 
enough) . 

The calculation of all ground states proceeds as fol- 
lows: Gi^i — Gi^i+i = for all feasible i holds. Start- 
ing with G,,j+2 = + 2)} {i = 1, . . . , L - 2) the 
complete set of ground states can be calculated recur- 
sively. The result is stored as a acyclic directed graph 
with Gij {1 < i < j < L) being the nodes and Gi^l 
the root. At each node, edges pointing to to the descen- 
dant sets Gi+ij , Gi+i^fc-i and Gk+i,i are stored instead 
of enumerating the states. Additionally, along with each 
node the ground-state energy _Eo(i, j) and the number of 
ground states di_j is kept. The degeneracy dij can be 
calculated recursively as well. 

The selection of a ground state is performed by a steep- 
est descent into the graph. Each ground state consists of 
the pairs encountered during the descent. At each node 
the steepest descent continues either into one descendant 


Gi+ij or into two descendant Gi+i^k-i,Gk+i,j, the al- 
ternative for proceeding being chosen randomly. The 
probability for each choice is proportional to the number 
of ground states found in the corresponding branch(s). 
For that purpose the degeneracy values di^ are used. 
It means that a path which contains twice the number 
of ground states of another path is selected on average 
twice as often. Therefore, it is guaranteed that each sin- 
gle ground state contributes with the same weight and a 
statistical correct T = average is obtained. 

For each sequence length the calculations were per- 
formed for 8000 independent realizations of the disor- 
dered, except for L = 2000 where only 1800 random se- 
quences were generated. For each realization 100 ground 
states were selected randomly and stored for further eval- 
uation. It was tested that by increasing this number the 
results do not change significantly. 

The resulting values a'^iq) are shown in Fig. |l| using a 
double logarithmic scale. Clearly, the function converges 
towards zero, thus P{q) ^ 5{q) for L — > oo, a similar 
result was found ^ for the model presented in For 
small sizes, this convergence is much smaller due to finite- 
size effects. This may be the reason that in [|] no decision 
about the behavior of the width of P{q) could be taken. 

The author thanks A. Pagnani and F. Ricci-Tersenghi 
for interesting discussions on the subject. 
A.K. Hartmann 

Institut fiir theoretische Physik, Gottingen, Germany 
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FIG. 1. Variance u^{q) of the distribution of overlaps as a 
function of system size L. The line represents the function 
O.TIL^""''. Please note the double logarithmic scale. 
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